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ANALYTICAL  MODEL  FOR  PREDICTING  THE  DYNAMIC  BEHAVIOR  OF  CAVITIES 


The  purpose  of  the  present  study  is  to  develop  an  analytical  model  for  predicting 
the  dynamic  behavior  of  cavities.  The  development  of  this  model  is  in  conjunction 
with  an  experimental  program  being  conducted  at  the  Applied  Research  Laboratory 
(ARL) , Pennsylvania  State  University. 


The  work  is  being  performed  under  funding  fran  Dr.  T.  E.  Peirce,  Naval  Sea  Systems 
Ccmnarxi  (SEA-03512)  (SR1230102,  element  6U52N). 
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INTRODUCTION 

The  objective  of  this  study  is  to  develop  an  analytical  model  for  predicting 
the  dynamic  behavior  of  cavities  formed  by  high-speed  underwater  vehicles  (vapor 
cavities)  and  cavities  formed  during  water  entry  and  underwater  launchings  (air/ 
vapor  cavities) . The  development  of  this  model  is  in  conjunction  with  an 
experimental  program  being  conducted  at  the  Applied  Research  Laboratory  (ARL) , 
Pennsylvania  State  University.  Recent  experimental  work  investigating  cavities 
formed  by  cones  in  a water  tunnel  at  ARL  is  reported  in  reference  (1) . The 
use  of  a water  tunnel  is  advantageous  in  that  the  mass  attrition  rates  can  be 
determined  for  different  size  cavities  by  measuring  the  air  flow  rate  feeding 
the  cavity. 


BACKGROUND 

The  need  to  predict  the  trajectory  of  cavity-running  underwater  vehicles 
has  been  given  emphasis  by  recent  development  effort.  To  successfully  predict 
the  trajectory  of  a vehicle  riding  in  a cavity  requires  a knowledge  of  the  forces 
acting  on  the  vehicle  which  are  influenced  by  the  presence  of  a cavity.  Water 
entry  and  underwater  launching  are  areas  of  weapon  delivery  where  the  presence  of 
a cavity  has  a significant  effect  on  the  trajectory  of  a vehicle. 

High-speed  underwater  vehicles  produce  a trailing  vapor  cavity  when  the 
minimum  local  pressure  on  the  vehicle  is  less  than  the  water  vapor  pressure.  At 
water  entry,  high-speed  vehicles  produce  an  open  air  filled  cavity.  Upon  cavity 
closure,  the  cavity  often  continues  to  expand  and  then  collapse  to  form  one  or  more 
stationary  cavities  and  a moving  trailing  cavity  which  behaves  in  the  same  manner 
as  a vapor  cavity.  In  fact,  this  trailing  cavity  can  beccme  a vapor  cavity  after 
the  entrapped  air  has  been  removed  by  attrition. 

The  pressure  within  these  cavities  and  their  size  and  shape  have  a significant 
effect  on  the  performance  and  trajectory  of  the  vehicle.  The  most  pronounced 
effect  is  on  the  drag.  When  a vehicle  is  trailed  by  a cavity,  the  drag  on  the 
vehicle  (other  than  friction  drag  which  can  be  significant  for  long  slender  bodies) 
is  primarily  due  to  the  dynamic  pressure  on  the  forebody  and  the  low  cavity  pressure 
on  the  afterbody.  There  is  no  pressure  recovery  on  the  afterbody  until  the  cavity 
closes  on  the  afterbody.  An  unstabilizing  lift  force  is  produced  on  the  forebody 
when  an  angle  of  attack  exists.  A lift  force  can  also  be  produced  on  the  afterbody 
when  it  contacts  the  cavity  wall.  Predicting  these  factors  can  best  be  accomplished 
by  first  predicting  the  behavior  and  configuration  of  the  cavity. 


Kim,  J.  H.  and  Holl,  J.  W. , "Water  Tunnel  Simulation  Study  of  the  later  Stages  of 
Water  Entry  of  Conical  Head  Bodies",  Applied  Research  laboratory,  Pennsylvania 
State  University,  Technical  Manorandum  75-177,  18  Jun  1975 
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ELLIPSOID  MODEL 


Munzer  and  Reichardt  (Ref.  (2))  and  Waid  (Ref.  (3))  described  the  cavity 
outline  as  a generalized  ellipsoid.  Munzer  and  Reichardt  derived  their  ellipsoid 
from  a theoretical  constant-pressure  surface  while  Waid  obtained  his  experimentally 
from  cavities  due  to  a disk.  A study  was  made  by  May  (Ref.  (4))  to  determine  how 
widely  these  two  representations  and  a true  ellipsoid  differed.  May  gave  a sketch, 
copied  in  Figure  1,  shewing  a comparison  of  the  generalized  ellipsoid  formulae 
using  an  axis  ratio  of  six.  The  outlines  are  indistinguishable  over  most  of  their 
length  and  differ  only  in  the  vicinity  of  the  nose. 


Reichardt  (Ref.  (5) ) derived  by  source-sink  methods  a theoretical  formula 
for  the  cavity  diameter,  D 


D 

d 


1/2 


(1) 


which  shews  the  diameter  to  be  a function  of  only  the  cavity-running  drag  coefficient  i 

and  the’  cavitation  number.  To  illustrate  why  the  form  of  equation  (1)  is  possible  i 

he  gave  a sketch,  copied  frem  reference  (4)  in  Figure  2,  which  shows  the  outlines  i 

of  the  cavities  formed  by  an  ellipsoidal  nose  and  a disk  with  the  same  drag  area, 

C jA.  The  cavities,  having  the  same  diameter,  differ  in  contour  only  in  the  ' 

vicinity  of  the  nose.  ; 

Since  an  ellipsoid,  for  which  potential  flew  solutions  are  available,  appears 
to  be  a reasonable  choice  for  a model  describing  the  cavity  shape,  the  following 
approach  was  taken.  It  is  assumed  that  the  pressure  drag  force  on  an  arbitrary 
nose  shape  which  forms  a cavity  is  equal  to  the  pressure  drag  force  on  the  fore- 
body of  an  ellipsoid  which  has  a minimum  local  pressure  equal  to  the  cavity 
pressure,  P . 


2 

Munzer,  H.  and  Reichardt,  H. , "Rotational  Synmetrical  Source-Sink  Bodies  with 
Predominantly  Constant  Pressure  Distributions",  Armament  Research  Establishment, 
Ministry  of  Supply,  Fort  Halstead,  Kent,  England,  Trans.  1/50,  1950 

^Waid,  R.  L. , "Cavity  Shapes  for  Circular  Disks  at  Angles  of  Attack",  California 
Institute  of  Technology  Hydrodynamics  Laboratory,  Report  E-73.4,  1957 

4 

May,  A.,  "Water  Entry  and  the  Cavity-Running  Behavior  of  Missiles",  Naval  Surface 
Weapons  Center,  White  Oak  Laboratory,  SEAHAC  TR  75-2,  1975 


^Reichardt,  H. , "The  Laws  of  Cavitation  Bubbles  at  Axially  Symmetrical  Bodies 
in  a Flew",  Ministry  of  Aircraft  Production  Volkennode,  MAP-VG,  Reports  and 


Translations  766,  ONR,  1946 
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The  local  pressure,  P,  on  the  surface  of  an  ovary  ellipsoid,  a > b 


(2) 


can  be  expressed  (Ref.  (6)) 


p - P 


-M^r) 


,i,2  a (4  - a ) 
(y  ) o o 


1 + (y*)2 


(3) 


where 


ao  = 2<V  ~ « 


^o  + 1 


f «oto  50  - 1 


- 1 


(4) 


£ 2 = 9 


o 2 2 

a - b 


(5) 


y'  " 

The  minimum  pressure,  P = Pc,  occurs  when  y'  = 0 


1 2 
P - P = i pU 
c °°  2 


1 - 


(dr) 


(6) 


(7) 


The  cavitation  number  is  defined 


P - P 


o = 


1 r,2 

2 PU 


, 2 

* - 
o' 


(3) 


^ Smith,  C.  W.  and  Guala,  J.  R.,  "Axisyrmetric  Ellipsoids  in  Ideal  Fluids", 
Naval  Surface  Weapons  Center,  White  Oak  Laboratory , NOI.TE  7 1- 1 00,  1 Ann  1973 
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and  the  local  pressure  coefficient  is  defined 


P - Pc 

Cp  = 1 2 

P 2 pU 


(1  + o)  (y')2 


Integrating  over  the  surface  of  the  forebody  of  an  ellipsoid,  o <_x<  a, 
the  drag  force  is  obtained 

b 

Drag  = 2tt  / (P  - Pc)ydy 


= 2tt  J [(P  - PJ  + (Pco  - Pc>]  Ydy 


1 TT2  ,2_, 

= 2 P°  1113  CD 


where  denotes  the  pressure  drag  coefficient  of  an  ellipsoid  based  on  the 

2 

maximum  frontal  area,  irb 


CD  - (1  + n)  - 1 


2 F ^ 

Cln_io.,_  _ 

L C - 1 


Equating  the  drag  on  the  nose  shape  forming  the  cavity  to  the  drag  on  the  ’cavity- 
ellipsoid' 

12  2 12  2 

Drag  = =■  rU  ib  C = ^ pu  nr  C,  (12) 

2 D 2 d 


The  follov;inq  expression  is  obtained  for  the  cavity  diameter 

V2 


. _^d_ 

Lcd  j 


where 


D = 2b,  cavity  diameter 
d - 2r , body  diameter 


* 


I 

f 

r 
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f 
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Table  1 

Comparison  of  Results 

Ellipsoid  Model  and  Reichardt's  Theoretical  Formula 


o 

Ellipsoid  Model, 
Equation  (15) 

Reichardt,  Equation 

.001 

1011. 

1015. 

.0021 

500. 

502. 

.0042 

250. 

251. 

.0069 

153.7 

154.4 

. 013C 

78.6 

79.0 

.0244 

44.2 

44.4 

.036 

30.2 

30.3 

.049 

22.1 

22.2 

.059 

18.4 

18.5 

.082 

13.5 

13.5 

.105 

10.5 

10.5 

.170 

6.55 

6.56 

.207 

5.39 

5.41 

.259 

4.32 

4.34 

.318 

3.52 

3.54 

.405 

2.77 

2.79 

.500 

2.24 

2.27 

.585 

1.91 

1.94 

.700 

1.60 

1.63 

.858 

1.30 

1.34 

.09 

1.02 

1.06 

i 

1 
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Rearranging  equations  (1)  and  (13)  as  follows,  for  equation  (1) 


a - .132o 


for  equation  (13) 


and  comparing  the  results,  it  is  seen  that  the  results,  as  shown  in  Table  1,  are 
nearly  identical,  especially  for  low  cavitation  numbers.  To  facilitate  the  cal- 
culations for  equation  (15)  the  axis  ratio,  a/b,  was  first  selected,  a was  then 
calculated  from  equations  (4)  and  (5),  the  cavitation  number,  o,  from  Bguation  (8) 
and  CD  from  equation  (11).  Since  the  agreement  between  equations  (14)  and  (15) 
is  so  close,  a more  direct  method  to  obtain  CL  can  be  formulated  by  equating 
equations  (14)  and  (15).  The  concept  of  the  ’cavity  ellipsoid',  however,  is 
retained  since  it  permits  a more  direct  treatment  for  the  dynamic  behavior  of 
the  cavity. 

Figures  3 through  6,  copied  from  reference  (4),  show  the  scaled  diameter  of 
cavities  plotted  from  data  obtained  by  Reichardt  (Ref.  (5)),  Waid  (Ref.  (3)), 
Rouse  and  McNcwn  (Ref.  (7))  and  Eisenberg-Pond  (Ref.  (8)).  Equation  (15), 
plotted  in  these  figures,  appears  to  be  a good  representation  of  the  data. 

Figures  7 through  9,  obtained  from  reference  (4),  show  data  for  the 
lengths  of  cavities.  The  cavity  half-lengths  were  measured  from  the  base  (or 
separation  point)  to  the  point  of  maximum  cavity  diameter.  The  ratios  of  the 
cavity  length  to  cavity  diameter,  D,  and  body  diameter,  d,  can  thus  be  expressed: 

L_L_x_xa  . 

D 2b  b a b (16) 


where  x/a  denotes  the  position  of  the  separation  point  on  the  cavitating  nose  within 
the  assumed  'cavity  ellipsoid'.  As  shown  in  Figure  2,  the  separation  point  is  not 
located  at  the  tip  (x/a  = 1)  of  the  'cavity  ellipsoid'.  In  correlating  the 
data  for  the  length  of  cavities  formed  by  a disk,  the  following  expression  was 
obtained : 


Rouse,  H.  and  McNcwn,  J.  S.,  "Cavitation  and  Pressure  Distribution;  Head  Forms 
at  Zero  Angle  of  Yaw",  State  University  of  Icwa,  Studies  in  Eng.  Bull.  32,  1948 


Eisenberg,  P.  and  Pond,  II.  L. , "Water  Tunnel  Investigations  of  Steady-State 
Cavities",  David  Taylor  Model  Basin  Report  660,  1948 
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(18) 


Using  equations  (18) , (13) , and  (5) , 


L 

D 


equations  (16)  and  (17)  can  be  expressed : 


2 

3 

2 

3 


a 


o 


- 1) 


a + (£  - 


1/2 


(19) 


_L = L 

DvC" 

d D 


(20) 


Equation  (20)  is  plotted  in  Figures  7 and  8 which  show  the  scaled  length 
of  cavities  as  a function  of  cavitation  number.  In  Figure  7 data  for  cavities 
due  to  a disk  are  plotted.  In  a similar  manner  data  for  various  nose  shapes 
are  plotted  in  Figure  8.  The  agreement  with  equation  (20)  is  shown  to  be 
excellent.  In  Figure  9 experimental  L/D  data  and  equation  (19)  are  plotted. 

The  agreement  with  the  experimental  data  is  shown  to  be  excellent  in  the 
region  for  a < 0.2.  In  the  region  o > 0.2,  equation  (19)  appears  to  be  an 
upper  bound  on  the  data  which  exhibits  considerable  scatter.  The  scatter  of 
the  data  in  Figure  9 may  be  due  to  the  del  xsndence  of  L/D  on  the  nose  shape 
(slope  at  the  separation  point  or  of  the  cavitatinq  nose) . 

Recent  experimental  L/D  data  obtained  at  ARL  (Ref.  (1))  for  18-  and  45-degree 
cones  with  cylindrical  afterbodies,  shown  in  Figure  10,  clearly  show  a strong 
dependence  on  the  nose  shape.  The  need  for  further  study  is  clearly  evident. 

The  dependence  of  the  cavity  L/D  on  the  nose  shape  suggests,  from  equation  (20), 
the  dependence  of  the  scaled  cavity  length,  and/or,  from  equation  (13) , the 
dependence  of  the  scaled  cavity  diameter  on  the  nose  shape.  The  aareement  with 
various  nose  shapes,  shown  in  Fiqure  8 for  the  scaled  length  of  cavities,  and  in 
Figures  3 and  5 for  the  scaled  diameters  of  cavities,  may  be  misleadina.  The 
experimental  data  in  these  fiqures  were  scaled  using  cavitating  drag  coefficients. 
Except  for  the  disk,  drag  data  are  sparse  in  the  region  for  ■ .2  and  drag 

coefficients  obtained  from  these  data  should  not  be  considered  reliable  enough 
to  establish  a firm  trend. 


DYNAMIC  MODEL 

The  velocity  potential  for  an  axisyrrmetric  ellipsoid  is  readily  obtained, 
as  shown  in  reference  (9),  using  dimensionless  ovary  semi-elliptic  coordinates 
(.i,  ').  The  transformation  equations  are: 


9 

Durand,  W.  F. , Aerodynamic  Theory,  Vol.  I,  Div.  C,  Chap  VI T,  Dover  Publications, 
Inc.,  New  York,  1963 
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where 


x = %/a2  - b2 


y = %/a2  - b2  Vl  - m2  %/t,2  - 1 


-1  < t)  < + 1 and  £ > 1 


the  coordinate  £,  which  can  be  written: 


(si+v 


, v > 0 


defines  the  surfaces  of  confocal  ellipses  with  the  foci  at 

x = ±ya2  - b2  , y = 0 (25) 

The  surface  of  the  'cavity  ellipsoid'  is  defined  by  the  coordinate  for  v = 0 


= f (v  = 0)  = 


L2  - b2 


and  the  Cartesian  coordinates  on  the  surface,  •'  , become: 

o 

x = a;.  (27) 

v = b - 2 (28) 

The  velocity  potential  function  for  an  axisyrmetric  ellipsoid  movinq  with  an 
axial  velocity,  U,  in  an  infinite  ideal  fluid  can  bo  expressed.  Ref.  (9). 


■l  1 Ar'  [5  ln  ’"H  “ '*] 


A.  = n V o"  - 1 «' 

For  an  oxnandinq  (or  contractinq)  ellipsoid  with  the  fluid  velocity  everywhere  normal 
to  the  confocal  olliosoids,  the  velocity  [xitontiul  function  can  lx?  expressed: 
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where 


and 


In 


S + 1 
£ - 1 


A2  = b b ? 


a ao 

ab  = ab  -*■  g-  = g-  , Constant 
o 


(31) 

(32) 

(33) 


The  velocity  potential  function  for  a moving  and  expanding  ellipsoid  can  new  be 
obtained  by  adding: 


<t>  = + 4>2 

The  fluid  velocity  corrponent  normal  to  the  oonfocal  ellipsoids  is: 


(34) 


and  the  tangential  fluid  velocity  carponent  is: 


(35) 


V = - 

y 


iF7 


{ 


l - y~ 

,2  2 

Z - y 


9<J> 

3y 


(36) 


On  the  surface  of  the  'cavity  ellipsoid'  the  fluid  velocity  ccnponents  are: 


(37) 


(38) 


To  determine  the  configuration  of  a dynamic  cavity,  consecutive  'cavity  disks' 
of  finite  thickness  can  be  tracted  in  inertial  coordinates  (Z,  R) , Figure  11.  Hue 
initial  properties  of  the  'cavity  disk'  (a  , b ) can  be  determined  using  the 
instantaneous  velocity  (UQ)  and  the  cavitypressure  (Pc  ) at  the  moment  (t  = tQ) 

the  fcavity  disk'  separates  from  the  cavitating  nose.  Tne  radius,  R,  and  position, 
Z,  of  the  'cavity  disk'  are: 
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= Y = b^l  - 


(39) 


-Zo  = x+/ 


Udx  = au  + / Udx 


/ 


(40) 


where  Z is  the  centroid  of  the  ellipsoid  which  contains  the  'cavity  disk'  formed 

at  t = 9. 

o 

Taking  the  time  derivative  of  equation  (40) : 


It  is  seen  that: 


and 


Z=x+U=ay+ay+U 


Z = 0 


x = ay  = -U 

when  the  cavity  is  not  expanding,  a = 0.  Thus, 


(41) 

(42) 

(43) 


• U 


and  for  an  expanding  cavity: 


(44) 


Z = ay 


(45) 


To  obtain  an  expression  for  b,  the  change  in  the  kinetic  energy  of  the  fluid 
flow  field  must  be  considered.  The  kinetic  energy  of  the  fluid  is: 


- 1 n f f cb  ^ 

- - 2 p J J ^ 


dS 


(46) 


K = €, 


where,  on  the  surface,  K - Kc 


dS  = 2ft (a2  - b2)  2 - y2)  (£  2 - l)du 

f o o 


30  _ ± 

9n  ~ JO  |T 


F~- 


- 1 3* 

▼ ~ T 


so  “ 11 


(47) 


(48) 
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The  change  in  the  kinetic  energy  due  to  the  expansion  of  the  cavity  can  be 
written: 

(P  - P ) ^ ~ (58) 

00  dt  dt 

where  P is  the  average  pressure  on  the  surface  of  the  'cavity  ellipsoid'.  Perform- 
ing the  differentiation  and  rearranging,  equation  (58)  becomes: 


P - P 


= K„ 


[b  b + \ b 2]  - | kxuz 


P 2 

Defining  the  cavity  pressure  when  o and  b = 0,  from  equation  (7) 


(59) 


P _ P 
co 


- s "‘[(^V)‘  - ■] 


(60) 


instead  of  using  the  average  pressure  on  the  surface  when  b and  b = 0,  from 
equation  ( 59 ) 


equation  (59)  becomes: 


P - P . 0 

— 1 = - I u2k. 

n 2 1 


P - P„ 


— -3L  = k9 

P 2 


[b  £ + | b- 2] 


(61) 


(62) 


Solving  for  b 


.2  p _ p 

b = -lL_+_ ^ 

2b  obK„ 


(63) 


an  expression  for  b can  be  obtained: 


t,=/£ 


dt 


(64) 


CAVITY  PRESSURE 

A cavity  pressure  model  for  predicting  the  cavity  pressure  and  mass  attrition 
rate  needs  to  be  developed.  The  model  will  need  to  incorporate  an  air  flew  model 
for  open  cavities  and  an  energy  model  for  expanding  cavities  which  are  often 
experienced  during  water  entry  after  cavity  closure.  An  attempt  to  correlate 
the  mass  attrition  data  (air  flow  coefficients)  contained  in  Reference  (1)  with 
one  correlation  for  18-  and  45-degree  cones  using  the  cavitating  drag  coefficient 
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as  a parameter  was  unsuccessful  due  to  the  lack  of  reliable  drag  coefficients  at 
high  cavitation  nunbers.  It  is  believed  that  the  mass  attrition  occurs  primarily 
at  the  base  (tail)  of  the  cavity  and  is  in  seme  manner  proportional  to  the  drag 
force  which  creates  the  cavity.  Correlation  of  mass  attrition  data  in  this  manner 
may  eliminate  the  need  for  a separate  correlation  for  each  drag  shape. 

RECCWMENDAT IONS  AND  CONCLUSIONS  i 

The  predictions  of  the  model  for  the  diameter  and  length  of  steady-state 
cavities  agrees  very  well  with  experimental  data  except  for  the  dependence 
of  the  cavity  L/D  on  nose  shape.  The  cavity  length,  as  shown  in  Figure  10 
and  also  noted  in  reference  (1) , was  found  to  be  strongly  influenced  by  the  body 
shape  at  high  cavitation  numbers.  The  drag  coefficient  is  a likely  parameter  for 
correlating  the  cavity  length  and  mass  attrition  rate.  Since  drag  data  at  high 
cavitation  numbers  is  sparse,  it  is  recommended  the  future  tests  conducted  at  ARL 
include  the  measurement  of  the  drag  force.  The  use  of  a water  tunnel  is  particularly 
appropriate  since  steady-state  conditions  can  be  maintained  during  the  measurement 
of  the  drag  force.  Future  tests  should  also  consider  additional  nose  shapes  with 
and  without  an  afterbody.  This  reconmendation  is  made  in  order  to  obtain  reliable 
correlations  which  hopefully  will  apply  to  any  drag  shape  (given  the  drag 
coefficient)  and  include  the  effect  of  cavity  closure  on  the  afterbody. 

Additional  work  on  the  analytical  model  is  required.  A model  for  predicting 
the  cavity  pressure  needs  to  be  developed  and  the  effect  of  the  free  surface  should 
be  included.  A digital  program  also  needs  to  be  written  and  the  predicted  behavior  ’ 

of  the  cavities  conpared  with  future  experimental  results  and  the  results  contained  | 

in  reference  (1).  The  experimental  results  contained  in  references  (10)  and  (11)  1 

for  the  cavity  shape  and  pressure  during  water  entry  should  also  be  compared  with  | 

the  predictions  of  the  analytical  model.  As  the  development  of  the  model  progresses,  , 

the  need  for  dynamic  water-entry  tests,  and  the  test  conditions,  can  be  better  j 

identified.  j 


°Aronson,  P.  M.,  "An  Experimental  Investigation  of  Cavity  Pressure  Scaling  and 
Drag  Coefficients  Encountered  in  the  Water  Entry  of  Vehic  les",  PhD  Thesis, 
University  of  Maryland,  1975 

11 

Abelson,  H.  I.,  "The  Behavior  of  the  Cavity  Formed  by  a Projectile  Entering  the 
Water  Vertically",  PhD  Thesis,  University  of  Maryland,  1969 


17 


TRUE  ELLIPSE 


NSWC/WOL/TR  76-114 


<D 


U> 


CO 


18 


FIG.  1 COMPARISON  OF  PROPOSED  CAVITY  OUTLINES,  REF  (4) 
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FIG.  2 SKETCH  OF  CAVITIES  DUE  TO  ELLIPSOID  AND  DISK  HAVING  THE 
SAME  DRAG  AREA  CdA.  AFTER  REICHARDT  REF  (5) 
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FIG.  3 CAVITY  DIAMETER  DATA  - REICHARDT  REF  (5) 
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FIG.  7 SCALED  LENGTH  OF  CAVITIES  DUE  TO  A DISK 
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FIG.  8 LENGTH  OF  CAVITIES  FOR  VARIOUS  NOSE  SHAPES 
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FIG.  9 EXPERIMENTAL  L/D  DATA 


26 


NSWC/WOI/TR  76-114 


DISTRIBUTION 

Copies 

Commander 

Naval  Air  Systems  Camand 
Washington,  D.  C.  20362 
Attn:  AIR-530133A 

AIR-53712A 
AIR-350D 
AIR-350 
AIR-350B 
AIR-030 

Commander 

Naval  Sea  Systems  Command 
Department  of  the  Navy 
Washington,  D.  C.  20362 
Attn:  SEA-09G32 

SEA-03  2 

SEA-0351 
SEA-03512 
SEA-0372 
SEA-054 
SEA-05121 
FHS-0647 
SEA-03B 

Office  of  Naval  Research 
Department  of  the  Navy 
Arlington , Virginia  22217 
Attn:  CKR- 4 38 

ONR-439 
ONR-463 

Director 

Defense  Research  and  Engineering 
The  Pentagon 
Washington,  D.  C.  20301 

Defense  Documentation  Center 
Cameron  Station 

Alexandria,  Virginia  22314  12 

Director  of  Research 

National  Aeronautics  and  Space  Administration 
600  Independence  Avenue  Southwest 
Washington , D.  C.  20546 


NSWC/WOL/TR  76-114 


Copies 


NASA  Scientific  and  Technical  Information  Facility 

Post  Office  Box  33 

College  Park,  Maryland  20740 

Superintendent 

U.S.  Naval  Postgraduate  School 
Monterey,  California  93940 
Attn;  Library 

Contending  Officer 
U.S.  Naval  Air  Development  Center 
Warminster,  Pennsylvania  18974 
Attn;  NADC  Library 

Director 

U.S.  Naval  Research  Laboratory 
Washington,  D.  C.  20390 
Attn:  Library 

Caimanding  Officer 

Naval  Ship  Research  and  Development  Center 
Bethesda,  Maryland  20084 
Attn;  Code  154 
Gode  1541 
Code  1556 
Code  5613 

Germander 

Naval  Undersea  Center 
San  Diego,  California  92132 
Attn:  Dr.  A.  Fabula 

Dr.  J.  Hoyt 

Ccrmanding  Officer 
Naval  Underwater  Systems  Center 
Newport,  Rhode  Island  02840 

Attn;  J.  F.  Brady 

W.  A.  McNally 

R.  H.  Nadolink 


Germander 

Naval  Weapons  Center 
auna  lake,  California 
Attn:  Library 


93555 


NSWC/WDL/TR  76-114 


i 


Harry  Diamond  Laboratories 
Washington,  D.  C.  20438 
f Attn:  Harry  Davis 

Library 

U.S.  Army  Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  Maryland  21105 
Attn:  Library 

National  Aeronautics  and  Space  Administration 
George  C.  Marshall  Space  Flight  Center 
Huntsville,  Alabama  35812 
Attn:  Library 

National  Aeronautics  and  Space  Administration 
Langley  Research  Center 
Langley  Station 
Hampton,  Virginia  23365 
Attn:  Library 

National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  Ohio  44135 
Attn:  Library 

Director 

Alden  Research  Laboratories 
Worcester  Polytechnic  Institute 
Holden,  Massachusetts  01520 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
Joins  Hopkins  Road 
Laurel,  Maryland  20810 
Attn:  Docum  ent  Librarian 

Applied  Research  Laboratory 
The  Pennsylvania  State  University 
Post  Office  Box  30 
State  College,  Pennsylvania  16801 
Attn:  Dr.  Blaine  R.  Parkin 

Dr.  J.  W.  Ho31 


AVCO  Fverett  Research  Laboratory' 
2385  Revere  Beach  Parlcv.ay 
Ew.rctt , Massachusetts  02149 
Attn:  Library 


Copies 


NSWC/VOI/TR  76-114 


Ba telle  Manorial  Institute 
505  King  Avenue 
Colimbus,  (Mo  43201 
Attn:  Library 

The  Boeing  Company  Aerospace  Library 
Post  Office  Box  3999 
Seattle,  Washington  98124 

Garfield  Thanas  Water  Tunnel  Library 

The  Pennsylvania  State  University 

Institute  for  Science  and  Engineering 

Applied  Research  Laboratory 

Post  Office  Box  30 

State  College,  Pennsylvania  16801 

General  Electric  Company 
Post  Office  Box  8555 
Philadelphia,  Pennsylvania  19101 
Attn:  MSD  Library 

Hydronautics,  Incorporated 
P indell  School  Road 
Laurel,  Maryland  20810 

Attn:  Mr.  Philip  Eisenburg 

Iowa  Institute  of  Hydraulic  Research 
State  University  of  Iowa 
Iowa  City,  Iowa  52240 
Attn:  Dr.  J.  F.  Kennedy 

Jet  Propulsion  laboratory 
4800  Oak  Grove  Drive 
Pasadena,  California  91103 
Attn : Library 

Lockheed  Missiles  and  Space  Company 
leper tment  57-24,  Building  150 
P.O.  Box  504 

Sunnyvale,  California  94C08 
Attn : Robert  Waid 

Ios  Alarms  Scientific  laboratory 
Post  Office  Bex  1G63 
Los  Alarms,  Few  Mexico  87544 
Attn:  Repor t Li bra ry 


4 


i 


MDW.NSWC<WVM05/1  (R««.  1-75) 


TO  AID  IN  UPDATING  THE  DISTRIBUTION  LIST 
FOR  NAVAL  SURFACE  WEAPONS  CENTER,  WHITE 
OAK  LABORATORY  TECHNICAL  REPORTS  PLEASE 
COMPLETE  THE  FORM  BELOW: 

TO  ALL  HOLDERS  OF  NSMC/WOL/TR  76-114 
by  Charles  W.  Smith  Code  WA-42 
DO  NOT  RETURN  THIS  FORM  IF  ALL  INFORMATION  IS  CURRENT 


A,  FACILITY  NAME  AND  AODRESS  (OLD)  (Show  Zip  Cede)  \ 

( 


j! 

NEW  ADDRESS  (Show  Zip  Code) 


B.  ATTENTION  LINE  ADDRESSES: 


C. 

Q1  REMOVE  THIS  FACILITY  FROM  THE  DISTRIBUTION  LIST  FOR  TECHNICAL  REPORTS  ON  THIS  SUBJECT . 


1 


DEPARTMENT  OF  THE  NAVY 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK.  SILVER  SPRING.  MD.  20910 


OFFICIAL  BUSINESS 
PENALTY  FOR  PRIVATE  USE.  S300 


POSTAGE  AND  FEES  PAID 
DEPARTMENT  OF  THE  NAVY 
DOD  316 


COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  SILVER  SPRING,  MARYLAND  20910 


ATTENTION:  CODE  WA-42 


■:kFHY*W  rtfei 


A 


